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ION TIME-OF-FLIGHT DETERMINATIONS' OF 
DOUBLY TO SINGLY IONIZED MERCURY ION RATIOS 
FROM A MERCURY ELECTRON BOMBARDMENT DISCHARGE 

J. M. Selleri, Jr., R. F. Kemp, and David F. Hall 
TRW Systems Group 
Redondo Beach, California 

I. INTRODUCTION 

The efficient operation of mercury electron bombardment ion thrusters 
for spacecraft propulsion requires a maximization of propellant utiliza- 
tion (i.e., the conversion of neutral mercury vapor to mercury ions) and 
an acceleration of the mercury ions to the optimum specific impulse for 
a given mission profile. Both of these goals may be adversely affected 
by the presence of multiply charged mercury ions in the thruster exhaust 
beam. For charge states higher than singly ionized mercury, the larger 
value of q/m for the particle (where q and m are particle charge and 
mass) leads to exhaust velocities which exceed the optimum specific 
impulse (assuming here that singly ionized mercury ions are accelerated 
to the optimum velocity) with a consequent loss of thrust efficiency. 

The presence of multiply charged ions in the exhaust beam also compli- 
cates the determination of propellant utilization when that utilization 
is based upon total exhaust ion current divided by the current equivalent 
of the injected neutral mercury vapor flow. Failure to correct the ion 
exhaust beam current flow for multiply charged ion species will lead to 
errors in the specification of propellant utilization for the thruster, 
with subsequent adverse effects on spacecraft mission performance. This 
program was undertaken to obtain data useful in establishing performance 
and test criteria for ion thrusters and in mission design. The resultant 
data have also provided a new approach to the diagnosis of electron 
bombardment discharge processes. 

The first experiments to determine the abundance of multiply charged 

ions in a thruster exhaust were conducted in 1962 by Milder^ and by Kemp, 

2 

Sellen, and Pawlik. Using a magnetic analyzer and hot wire current 
detecting probes. Milder obtained Hg ++ /Hg + ratios (hereafter denoted as 
++/+ ratios) as a function of the anode potential in the electron bom- 
bardment discharge. The onset of doubly charged ions was observed as 
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discharge potential was increased above ^20 volts, afid ++/+ ratios as 
large as 0.2 were obtained for anode potentials of 100 volts. The 
abundance of Hg ++ relative to Hg + for these experiments appeared to be 
larger than would occur if only single stage processes occurred (i.e., 
the production of Hg + and Hg ++ by a single ionizing collision of an 
electron) and suggested that "sequential processes" were at work (i.e., 
the production of Hg ++ by successive ionization encounters each of which 
resulted in the removal of a single electron). Milder also searched for 
Hg +++ , but was not able to detect this species, noting, however, that 
resolution of the magnetic analyzer was such that small quantities of the 
triply charged ion could be present without system detection. 

Kemp, Sellen, and Pawlik in their early experiments detected Hg ++ 
by two methods. The first method utilized a pulsed ion engine accelera- 
tion voltage with the ion current to a distant collector recording the 
arrival of, at first, the more rapidly moving doubly charged ions, and, 
after an appropriate delay, the arrival of the more slowly moving Hg + 
ions. The ++/+ ratio was approximately 0.06 which, for the discharge 
potential of 50 volts, was in comparatively good agreement with the 
results of Milder. The second experimental treatment of the ++/+ ratio 
utilized a steady-state ion exhaust beam flow into an ion velocity analyzer 
(operating as a shortened Wien velocity filter). Abundance ratios from 
this analyzer were in agreement with those obtained by the pulsed beam 
ion time-of-fl ight technique. When the discharge potential was increased 
to 100 volts, the pulsed beam tests revealed the presence of Hg +++ . 

3 

More recently, Masek performed exploratory measurements on a mercury 
ion engine exhaust with a variable field magnetic analyzer. Doubly 
charged ions of mercury were observed, but system resolution was not 
sufficient for other than order of magnitude measurements. The improve- 
ment of resolution in this experiment has led to the current findings by 

4 / 

Pawlik, Goldstein, Fitzgerald, and Adams. Their measurements (as had 
those of Milder) revealed an onset in Hg production for anode discharge 
potentials in excess of ^23 volts, and were interpreted as suggesting 
that Hg ++ is not produced, in the main, by single stage ionization pro- 
cesses from the Hg° ground state. Their experiments also revealed 
significant variations in the -H-/+ ratio as a function of position over 
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the face of the accelerator grid of the ion engine with predominant ++ 
production on the axis of the thruster and reduced ++ production near 
the boundaries of the engine. These measurements were also interpreted 
as evidence that sequential ionization to the doubly ionized state is 
the principal means for the production of this ion species. 

The experiments and analysis to be reported here will both continue 

and extend the previous investigations. The ion engine in use in the 

experiments is a 20-cm-diameter electromagnet electron bombardment ion 

4 

engine and is similar to that utilized by Pawlik et al . The experi- 
mental technique to determine ++/+ abundance ratios utilized ion time- 

of-flight, but was a variation of that technique used in the earlier 

2 

tests of Kemp, Seller, and Pawlik. There, the entire ion beam was pulsed 
on rapidly. In the present experiments, the ion flow from the ion engine 
is constant. A section of that ion flow enters an analyzer and is acted 
upon there by a transverse electric field which bends the ions. The 
rapid release of this deflection field allows the undisturbed ion flow to 
proceed directly into a Faraday cup with ion time-of-fl ight now occurring 
within the confines of the analyzer. This technique will be discussed 
further. In addition, an analytical treatment of the discharge action in 
producing various ion species has been carried out. That analysis is 
presented in the following section, and will be followed by a description 
of the time-of-fl ight technique, experimental results, and a discussion 
of those results. 
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II. ANALYSIS OF DISCHARGE PROCESSES 


In his treatment of Hg + and Hg ++ from electron bombardment discharges. 

Milder 1 identified three reactions as the basis of a simple analysis. 

4 

Pawlik et al broadened the number of reactive processes to be considered 
to include ionization of both atoms and ions from states other than the 
ground state. In a more recent work, Wu has numerically evaluated 
reaction rates for a series of possible ionization processes. In these 
calculations, it has been assumed that electrons possess a Maxwellian 
energy distribution. 

The following analysis differs somewhat in concept from previous 
approaches. Its basis is a consideration of a group of neutral mercury 
atoms, of number n Q , as that group moves through the electron bombard- 
ment discharge. In an actual ion engine, neutrals moving through the 
discharge undergo various changes of state in both ionization and 
excitation. The particles so formed may move to the exit plane of the 
engine, pass through the openings in the screen grid and undergo accel- 
eration. They may also intercept warm material surfaces within the 
discharge chamber. It is assumed that interception with such surfaces 
returns the mercury particle to the ground state of the neutral atom, 
whereupon it re-evaporates and re-enters the discharge. At some point 
in the particle "history" in the discharge, the particle trajectory will 
result in passage through the accelerator grid. The analysis will focus 
attention on that portion of the particle path between its last inter- 
ception with a material surface and its passage through the accelerator 
grid. This geometry is shown in Figure 1. 

Considering a group of n Q atoms and denoting t = 0 as the instant 
of the last contact with a material boundary before passage through the 
accelerator, the following equations may be written: 

dn 

Ht~ = ' F 01 n e"o - F 02 n e"o <’> 

dn + 

3t"~ = F 0l"e"o ' F 12 n e n + < 2 > 
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Figure 1. Simplified Discharge Configuration Considered in the Analysis. 

At t = 0 the neutral mercury atom completes its final encounter 
with a material boundary before passage through the accelerator 
grid at t = T. Excitation and ionization of the particle occurs 
along the trajectory, resulting in various abundance levels of 
Hg° , Hg + , and Hg ++ . 


= F 02 n e n o + F 1 2 n e n + 


The quantities n g , F Q1 , F Q2> and F 12 are defined below. Equations (1), 
(2), and (3) have considered three possible reactive processes: 

e + Hg° — ►Hg + + e + e (4) 


e + Hg° — ►Hg ++ + e + e + e (5) 

and . ++ 

e + Hg Hg + e + e (6) 


In the description here no assumption will be made concerning the 
excitation states of Hg°, Hg + , or Hg . Thus, metastable levels of 
Hg° are included in the process given by Equations (4) and (5). 
Similarly, the ions on the left hand side in Equation (6) may be in 
either the ground state or excited states. 
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In writing Equations (1), (2), (3) as the sole Reactive processes 

for treatment, the analysis has not considered electron- ion recombination 

reactions or charge exchange processes (while charge exchange between Hg° 

and Hg + does not result in a change in specie abundance, charge exchange 

between Hg ++ and Hg° which leads to 2 Hg + would alter such abundances). 

Also neglected here are the third and higher steps in the sequential 

-3 

ionization chain. The term n g is electron density (cm ) and it is 

important to note that, in an actual discharge, major changes will occur 

in this variable as particles move along their paths (increasing t). In 

addition to changes in the magnitude of n g , the energy distribution of 

electrons also changes from point to point in the discharge. The energy 

distribution function for electrons will be denoted as f where 

e 

/f e (E) dE = 1 (7) 

o 

The final terms to be defined are F Q1 , F Q2 , and F ]2 . In general, 

/ co 

f e< E > v e 5 ij< E > dE < 8 ) 

0 

is a net transfer probability (in cm 3 /sec) from the i th to the j th state 

3 

of ionization. When multiplied by n Q (electrons/cm ) n F. . has the 
- 1 ® ® i J 

dimensions of t . Note that the definition of F. . utilizes 5 - • as a 

t h ^ t h ^ ^ 

cross section for the transition from the i to the j ionization 

state. To make Equation ( 8 ) completely general a., includes transitions 

th ‘ J ■ 

from all possible excitation levels in the i state (thus Oq-j includes 

ionization from metastable neutral states). The use of the a-, notation 

* J 

is meant to indicate an appropriate weighting of the cross sections to 
be utilized over the distributions of n Q and n + in their initial states 
of the reaction. 

Returning to Equations (1), (2), and (3) and with the specification 
of Equation ( 8 ), it is possible in principle, to evaluate specie abundance 
from the original group of n Q atoms through integration along the particle 
path through the discharge (note that an integral in t may be transformed 
into an integral in s, distance along the particle trajectory, through use 
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of the relation ds = v dt). Finally, to evaluate the total ++/+ ratio 

r 

in the emerging ion beam, a weighted average would be required for each 
exit hole in the accelerator grid, with appropriate weight to all of the 
possible particle trajectories which ended in passage through the hole, 
followed by a summation, with appropriate weighting in current/hole, over 
the totality of holes in the exit grid. 

While the procedure outlined would result in an evaluated ++/+ ratio 
for the exhaust beam, the estimated extent of required numerical calcula- 
tions are considerable and beyond the intent of the present analysis. 
Variations in f and n occur along the particle path, and considerable 

G 6 

uncertainty exists over certain elements of the terms in 5... From the 

"I J 

measurements of Masek, electrons in the discharge consist of both "high 
energy primary" electrons and a Maxwellian group of "secondary electrons" 
with radically differing ionization probabilities because of the strong 
energy dependence in those terms. Moreover, from the recent measurements 
of Wilbur/ plasma potential variations in the discharge regions are of a 
complicated nature with accompanying complications in electron density 
and temperature. The potential variations in the plasma discharge impact 
upon integrals in time along the particle paths because of velocity 
changes of charged particles in the discharge electric fields (dt and ds 
are linked by v p , particle velocity). Of even further concern but not, 
as yet, treated, are particle orbit changes due to the turbulence of the 
electric fields in the plasma discharge. 

In view of the many real difficulties in assessing the ++/+ ratio 
by a rigorous evaluation of Equations (1 through 3) with appropriate 
weighting, an analytical approach would appear to be unwarranted. However, 
some qualitative insight into discharge behavior might be gained by a 
simple integration of these equations. The results could then be compared 
to the results of experiments to note if discharge behavior can be pre- 
dicted on, at least, a parametric basis. When these integrations are 
carried out the results are: 


n 0 (T) = 


n 0 (0)e 


‘ n e^ F 01 


f o 2 ) t 


(9) 
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( 10 ) 


n o^°) F 01 


-r, Q (Fm + F n5 )T " n o^i9" F 

„ /t\ . u ui / „ e v 01 02' a e iz 

n + ( T ) = r rr — +r — r * e - e 

+ f 12 - (F 01 + f 02 ) 


rT x n o^°^ F 02 (, " n e (F 01 + F 02^ T ) 

n++(T) ' F oT rrT o? 0 ' e > 


n o (0)F 01 F 12 

(l - e"^ 01 + F ° 2)T 

).( 


F 12 ‘ ^01 + F 02^ 

F 01 + F 02 


F 12 


(ID 


These straightforward results are, of course, only obtained in the 
case where n g and are invariant along the path of integration. With 
this as the principal simplifying assumption, estimates of the ++/+ ratio 
can be obtained for comparison to experimental findings. 

Note that the abundances are given in terms of the particle time of 
flight, T, from the last encounter with a material boundary and until 
passage through the acceleration grid. For the exhaust beam of an ion 
engine, the value of T will be given approximately by L/v p where L is 
the final drift path length in the discharge region and v p is particle 
drift speed averaged over the passage through the discharge. Consider- 
ing a drift distance of ~10 cm, a particle having 0.05 ev (~500°K) 

4 

thermal energy would travel at -vlO cm/ sec and escape through the accel- 
-4 -3 

erator grid in 10 < T < 10 sec. To evaluate the ++/+ ratio. Equation 

(11) is divided by Equation (10) to yield 




" n e (F 01 +F 02^) 


( 12 ) 




n e^ F 01 + F 02^ T " n e F 12 T 
e - e 


Equation (12) is not easily visualized in general in terms of ion 
abundance. However, some simplification may be made in the context of 
conventional ion engine operating conditions. Ion engines, if they are 
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to be utilized as thrusting devices, must have a high propellant util i 
zation efficiency. Since the emerging neutral fraction of the initial 


n e ( Fq-j + Fn<? ^ 


02 - 


neutral input is given by e 

requires that the discharge conditions be such as to yield |n g (F 
of at least 2. 


, a high propellant utilization 

+ F 02> T 


01 


An alternative statement of this condition is that the discharge 
must be "thick" to the removal of mercury from its neutral state. A 
second implied condition in efficient ion engine operation is that the 
discharge must be "thin" with regard to delivery of the particles into 
the doubly ionized state. These conditions are, in fact, obtained in 
the present development of ion engines, and, with these known indications 
on the magnitude of particular ionization levels. Equation (12) may be 
simpl ified to 


n F, J - e n - 6 n(e - 1 ) 
++/+ = - 5 — !± 

n - " e F 12 T 


(13) 


n 0 (°) " n 0 ( T ) 

where q is the fractional propellant utilization efficiency ( fTTOl 

P p 0 

r ip r np 

e = ■= — p and 0 = p^. In the limit of propellant utilization near 
h 01 h 02 h 01 

unity this expression becomes 


++/+ = 


n e F l 2 T 


e + 0 - £0 

n e F l 2 T 


(14) 


This is the expression which shall be used in an examination of experi- 
mental results. It should be noted that both 0 and e are energy 
dependent terms. Increases in anode discharge potential will presumably 
lead to a more energetic electron distribution for which both F Q2 and F 12 
will increase relative to Fq-j thus increasing both e and 0. It should 
also be noted that the only plasma density dependent term in the ++/+ 
ratio appears in n e F-| 2 T. Thus, variations in the discharge which result 
in increases in plasma density and which do not result in alterations 
in the electron energy distribution should affect the ++/+ ratio through 
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the term, and experiments revealing such density dependent 

behavior would be evidence that sequential ionization does occur in 
the discharge (non-zero F-| 2 ). 

It is also of interest to examine Equation (14) in the limit that 
certain transitional probabilities vanish. For example, if F-^ is set 
equal to zero, then doubly charged ions appear only through the 0— ~2 
transition. Examining (14) in this limit 


++/+ 





(15) 


This, it may be noted recovers the situation originally treated by 
Milder [Equation (9), Reference 1]. 

To illustrate the agreement between exact and approximate forms of 
the ++/+ ratio. Equation (12) has been evaluated as a function of n e /n e , 
electron density normalized to the value for which n g F 01 T ■ 2 - 
F q2 = 0.01 F Q i and F ]2 = 0.1 F Ql . These calculations are illustrated in 
Figure 2. As n g /n e — 0, the ++/+ ratio goes to 0.01 = f q 2 /Fqi . Thus 
++/+ — F q2 /Fq-| as n g — 0, irrespective of the magnitude of F 12 , and, also 
from Equation (15), as F^ — 0, irrespective of the magnitude of n g . Note 
that Equation (13) cannot be used for expressions of ++/+ in the limit 
of n g —0, since this approximate form of Equation (12) has a required 
condition of high propellant utilization which fails in the limit of 
vanishing electron density. For conditions in which propellant utiliza- 
tion is high. Equation (14) is in good agreement with the exact expression 
[Equation (12)] as may be noted in the curves of Figure 2. 

Another condition to examine is the limit of Equation (14) in the 
limit F Q2 — This condl ' t '' on is a comparatively likely one for discharge 
anode potentials less than <35 volts. Setting F Q2 to zero yields 


++ 

+ 


lim 

F 02+0 


n e F 12 T ~ 

1 " n e F 12 


e 

T 


(16) 
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Figure 2. The ++/+ Abundance Ratio versus Normalized Electron Density. 

Normalization is to n g , that density for which fi e F 01 T - 2 ' 
F 12^ F 01 = snd F 02 / F 01 - 0.01. 


where e = 




now becomes 



For n e F 12 T « 1, (a discharge "thin" with respect to the second step in 
the sequential ionization), then (16) may be approximated by 


++/+ = n e F-| 2 T 0 + n e F i2 T )^ " 


n e F 01 T ' 


07) 


Since the discharge in the engine is "thick" for the first ionization for 
cases of practical interest, the term n e FQ-|T is "large" and is essentially 
invariant for changes in, for example, discharge anode potential. 

Equation (17) may then be further approximated to 


and 


++ /+ = n e F i2 T ( 1 + n e F 12 T ) 


(18) 


3 (++/+) 
3v anode 


3 ( n e F 12^ 
3v anode 


0 + 2 n e F 12 T) 


(19) 
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Thus, if the anode potential can be varied without variations being 
introduced in plasma density, then electron energy distributions will 
be altered, and F-j g will be altered. Since F-^ rises with increasing 

electron energy, , should be expected to be > 0 and 

u anoae dV anode 

is expected to be positive. 

In conclusion to this present discussion, it should be emphasized 
again that the purpose of the analysis is to gain insight into ion engine 
behavior in at least its parametric dependences. It is recognized that 
the ion engine and its discharge is, in the overall sense, a complicated 
process. Analyses of ion engine propellant utilization, such as those 
by Kaufman and Kaufman and Cohen indicate that neutral propellant losses 
at high discharge loss conditions are approximately constant, irrespective 
of the total neutral injection flow rate. From this evidence, there are 
processes present in the engine discharge which are not represented in 
the original discharge Equations (1), (2), and (3), or if present, are 
subject to subtle changes as the neutral injection rate is varied. The 
analysis in this present treatment, then, operates with clearly acknowledged 
limitations, but with the possibility of some practical application. 
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III. ION TIME-OF-FLIGHT (TOF) TECHNIQUE 


A. OPERATIONAL PRINCIPLES AND PROBLEMS IN ION TOF ANALYZERS 

The principles of ion time-of-fl ight as a means of identifying ion 
charge species in the exhaust beam of an ion engine, were first used by 

p 

Kemp, Sellen, and Pawlik. In those experiments the acceleration voltage 

of the ion engine was pulsed on rapidly with all other thruster voltages 

and currents on and in steady state. Multiply charged species acquire 

1/2 • 

an exhaust velocity proportional to q where q is ion charge with a 
resultant separation of species in the initial TOF down the testing 
chamber. The ion current arriving at a distant collector has the general 
characteristics shown in Figure 3 (from Reference 2) (which has assumed 
only ++ and + species from the ion engine discharge). Separation of 
species, even if Hg +++ is present, is extremely clean provided that 
variations in ion time-of-flight from one charge species to another are 
large compared to the acceleration voltage switch-on time ('vl x 10 ^ 
seconds) . 

While the pulsed beam technique utilized in these previous experi- 
ments operated successfully, it is believed some concern could be raised 
over possible changes in the discharge operating conditions for the 
manner in which the pulsed technique was applied. Although there is no 
evidence for ion engine alterations, changes in the discharge might, in 
principle, occur depending on whether ions arriving at the screen grid 
were subject to extraction and acceleration or were allowed to be inter- 
cepted by and terminated at that electrode. Accordingly, for this series 
of experiments a variation in the ion TOF technique was utilized. In 
this variation the ion engine discharge and ion flow operate continuously 
and pulsed voltages are applied to a small portion of the ion exhaust 
beam approximately 1 meter downstream from the exit plane. The ion 
current entering a small aperture in the analyzer is subjected to 
electric fields which prevent the ions from reaching the analyzer collector. 
When these fields are removed, the now unperturbed ion flow moves through 
the analyzer, with, again an identifiable separation of specie currents 
during the initial ion TOF in the analyzer. The requirement for a clean 
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Figure 3. Representative Ion Current Signal to Collector Ring 1 in Pulsed 
Beam Tests After Acceleration Voltage is Pulsed on (t = 0). 

Beam Current component of Hg arrives at ^29 ysec after t = 0 
while beam current component of Hg + arrives at %41 ysec, with 
ion time-of-fl ight proportional to q -1 /2 as predicted. Ratio 
of J++/J+ is ^0.12, and relative abundance of ++/+ species is 
^0.06. No evidence of Hg +++ for discharge conditions utilized. 

2 

From Kemp, Sellen, and Pawlik 


separation of species is, in a manner similar to the previous TOF technique, 
that the switch-off time for the electric fields preventing ion passage in 
the analyzer be small compared to variations of ion TOF in the analyzer. 

For the principle of ion TOF considered thus far, two major possi- 
bilities exist in the "perturbation" electric fields in the analyzer. 

In one configuration these electric fields are "longitudinal" and are 
oppositely directed to the ion flow. Figure 4a illustrates this configu- 
ration. With sufficient retarding potentials, the previously accelerated 
beam ions will come to rest and then be reflected backward toward the 
engine. Removal of the retardation voltage allows the ions to continue 
their flow into the analyzer and their eventual detection at the analyzer 
collector. This method presents several problems. The first of these is 
that the potential applied to analyzer electrodes must be such as to bring 
the accelerated ions to rest. For typical ion engine operating conditions, 
the ion acceleration voltage is ^ 2000 volts and these potentials must also 
appear on the repeller grid in the analyzer. Such positive potentials. 
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a. LONGITUDINAL FIELD ION TOF ANALYZER 



INLET GRID WITH 



FARADAY CUP 
DETECTOR. 


b. TRANSVERSE FIELD ION TOF ANALYZER , ON AND ELECTRON 

TRAJECTORIES WHEN 
DEFLECTION VOLTAGE 

ELECTRON AND ION REMOVED 



Figure 4. Operational Principles of Longitudinal and Transverse Field 
Ion TOF Analyzers 

while repelling ions, attract electrons from the neutralized thrust beam 
and create electron drainage conditions which are highly disruptive to 
the beam neutralization. Accordingly, and to prevent electron drainage, 
an inlet grid on the analyzer must be used which will prevent electron 
entry into the analyzer, even in the presence of the comparatively strong 
attractive electric fields (for electrons) which exist due to the positive 
potentials on the ion repel ler. The introduction of grids into the ion 
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flow, however, is, in general, a perturbing force. Wi th the removal of 
the ion voltage on the repeller electrode, the ions may now proceed down 
the analyzer. However, electrons are still denied entry into the analyzer 
because of the inlet grid, and ions in the analyzer must proceed without 
the presence of neutralizing electrons. Since space charge forces in ion 
flows 6f this density are appreciable even for considerable downstream 
separations from the engine, ion beam expansion will occur in the analyzer 
with effects upon analyzer results which are difficult to assess. 

The alternative arrangement of electric fields in an ion TOF analyzer 
is a "transverse field" analyzer (Figure 4b). In this configuration the 
ion flow enters an inlet channel which imposes an electric field perpendi- 
cular to the ion flow direction. The ion is deflected in this field and is 
either collected on one of the inlet plates or is allowed to intercept the 
wall of the analyzer. The removal of the deflection voltage allows the 
ions to proceed with unperturbed trajectories to the collector, with, again, 
charge specie currents separating on the initial TOF in the analyzer. 

Several advantages present themselves in this configuration. The first of 
these is that deflection of the ions can be accomplished, in general, with 
potentials that are small compared to those required for ion reflection. 

A second advantage is that there is no requirement for positive potentials 
to be applied to the deflection plates. Ion deflection can be achieved 
with one plate at zero potential (neutralized beam potential*) and the 
remaining plate at a negative voltage. Ions are attracted to the negative 
plate and electrons are repelled by it. Plasma flow in the presence of 
such negatively biased electrodes remains stable (whereas the flow becomes 
highly disturbed and unstable for positive voltages: for discussion of 

these effects see, for example. Reference 10). With the removal of the 
negative deflection voltage, the unperturbed plasma flow begins its passage 
through the analyzer. Because grids are not required at the analyzer input, 
electrons may now accompany the ions so that ion motion in the analyzer is 
unperturbed by unneutralized space charge forces. Separation of electrons 
from ions in the flow does occur eventually but only as a result of bias 

*It is important to maintain all surfaces near the entrance aperture at 
this potential to minimize electrostatic "lens" effects which cause 
beam dispersion within the drift tube. 
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grids over the surface of the Faraday cup collector of the analyzer. 

From these several configurational advantages, a "deflection" type 
analyzer was chosen, rather than the "reflection" type. 

B. DEFLECTION VOLTAGE REQUIREMENTS IN A TRANSVERSE FIELD ION TOF ANALYZER 

Figure 4b has illustrated a deflection channel of width w and length 
X with a deflection voltage, V d , imposed across it. For an ion of mass, 
m, charge, q, and acceleration velocity, v + , the action of the deflecting 
electric field is to give the ion a transverse velocity, v t , in its 
passage through the channel. Upon emerging from the deflection plates 


tr 


iv 

m v + w 


( 20 ) 


and, since 


2 V. 


1/2 


»♦ ■ <TT> 


where V Q is the ion acceleration voltage, the angular deflection 6 upon 
emerging from the deflection plates is 


e 



( 21 ) 


A minimum requirement In the TOF analyzer is that an ion be deflected 
out of the Faraday cup at the end of the analyzer. For a cup of diameter 
D and with a deflector-to-collector separation of A, the minimum permissible 
value of e is 

« _ D _ X V d,min 

0 min 21 “ 2(TV^ (22) 

which yields 

w D V 

V d,min ” A A (23) 


Note that this is independent of ion charge, q, which is a desirable 
feature. 
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Now, as an example of the voltages required for- a practical case, 
consider w = 0.01m, D = 0.05m, X = 0.03m, and A = 0.6m with V Q = 2000 
volts. Then V d - n = 55 volts which is « V Q . Thus, only comparatively 
small deflection voltages are necessary for realistic experimental con- 
ditions in analyzer construction and ion energy. However, in practice 
it is desirable that ions be bent through angles which are larger than 
the minimum angle and, hence, result in interception on the analyzer 
support walls well in advance of the Faraday cup position to accommodate 
any beam dispersion. Also the above equations do not treat the effects 
of space charge on the required deflection voltage (discussed below). 

For both of these reasons additional deflection voltage will be required. 
Again, it is of interest to note that in the equations of motion, the 
ion charge q does not appear in the final ion trajectory. Thus, the 
transverse electrostatic deflection acts upon Hg + and Hg ++ in an identical 
manner, without the establishment of varying deflection requirements for 
the varying charge species. 


The problems of "space charge loading" are introduced when realistic 
values of ion density are treated. [Equations (20) through (23) represent 
the flow of a single ion]. In order for a transverse electric field to 
act upon all ions, it is necessary for V d to exceed the potential buildup 
in an unneutralized slab of ions of charge density p + and thickness w. 
Using MKS units, this requires 



(24) 


I C O 

For example, an ion flow density of 10 ions/m (and assuming singly 
charged ions) and for w = 0.01 meter, then 

V j • = 1.6(10)" 19 (10) 15 (10)" 4 (2 tt)( 9)(10) 9 = 900 volts, 

Q 9 III I M 


a voltage considerably in excess of that given from Equation (23) 

Since the deflection requirement now includes the effects of ion 
space charge, it is of interest to evaluate this quantity. In the plane 
immediately following ion acceleration and for 1 ampere ion flow moving 
at 4 x 10 4 m/sec and a beam of diameter 0.2 meter, the ion density is 
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= 5 x 10 15 ions/m 3 , which is five times larger than the example given 
above. The transverse deflection, then, is not appropriate if the 
analyzer is to be placed immediately downstream from the ion engine 
exhaust plane. However, because of angular divergence of the ion beam, 
the density at a point 1 meter downstream from the ion engine is reduced 
by an order of magnitude from its initial density. In these experiments 
the analyzer inlet was at 1 meter, where the requirement on due to 
space charge reduced to a value of ~ 450 volts for even the largest 
accelerated ion flows. In practice, values of employed ranged to 
^600 volts with the upper part of this range used for the highest levels 
of ion density. 

A final area of concern is in requirements on switch-off time for 
the deflection voltage and in the separation ability of the TOF analyzer 
in the identification of charge species. Since the deflection channel 
has, of necessity, a non-zero length, X, ions passing through it have a 

\ 4 

dwell time of For X = 0.03 meter and v + = 4 x 10 m/sec, this dwell 

v + 

time is approximately 1 microsecond. 

Thus, removal of the deflection voltage instantaneously would lead 
to a portion of the flow <vl microsecond in duration which would possess 
varying "histories" of transverse acceleration within it. On the other 
hand, if the deflection voltage switch-off occurs over a period long 
compared to 1 microsecond, then the Ion behavior downstream of'the Inlet 
channel will possess time variations mainly Influenced by the switch-off. 

An acceptable switch-off criterion, then, is to make the switch-off time 
less than ion dwell time in the deflection channel and to require, finally, 
that ion dwell time in the inlet be small compared to variations in the 
ion TOF in the analyzer due to variations in ion charge and ion exhaust 
velocity. For a deflector-to-col lector interspace of 0.6 meter and 
v + = 4 x 10^ m/sec, the x + (for Hg + ) is “15 microseconds. The t ++ (for 
Hg ++ ) is reduced to =10.5 microseconds (since v ++ = n/2v + ) and t +++ is 
=8.5 microseconds. The Ion dwell time in the inlet channel of M 
microsecond is small compared to A(t + - t ++ ) of 5 microseconds, but is of 
comparable magnitude to A(t ++ - x +++ ) of 2 microseconds. Therefore, 
separation of Hg ++ from Hg + should be well defined while separation of 
Hg ++ and Hg +++ would be marginal in the present system. However, from 
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the earlier measurements of Reference 2, the "appearance" potential in 
anode discharge voltage for Hg production is considerably above the 
parameter range for present ion engine operation. If, however, conditions 
should alter somehow so as to require these higher anode voltages and if 
the Hg +++ presence were of concern, the ion TOF configuration could be 
suitably altered to determine Hg +++ by increasing A or by decreasing V Q , 
the ion engine acceleration voltage. For the purposes of the present 
experiment, good separation of Hg ++ and Hg + was achieved with an analyzer 
drift space of 0.6 meter and an inlet channel length of 0.03 meter, 
deflection voltages of less than or equal to 600 volts, and switch-off 
times of less than 1 microsecond. 

C. ION TOF ANALYZER CONSTRUCTION AND OPERATION 

The values of A, A, w, and D have been previously specified, while 
Figure 4b has illustrated these quantities relative to the ion flow. 
Figures 5 and 6 show portions of the completed probe. In Figure 7 the 
arrangement and range of angular positions of the probe relative to 
the ion engine may be seen. 

Figure 8 illustrates the timing sequence for the switch-off of the 
acceleration voltage and Figure 9 shows oscilloscope traces with indicated 
arrivals and abundances of the Hg ++ species and the Hg + species. After 
these initial experiments, the TOF analyzer was then applied to the series 
of experiments described in the following section. 
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Figure 5. Deflection Plate Assembly and Faraday-Cup Assembly of the 

TOF Probe 



Figure 6. Completely 
Assembled Time-of-Flight 
Probe for Measurement of 
Hg + and Hg ++ Ion Abun- 
dances. The figure shows 
the 1 x 3 cm ion entrance 
aperture and clamps which 
affix the probe to the 
support structure. 
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Figure 7. Arrangement of Ion Engine and TOF Probe in Vacuum Tank. 

Engine projects beam vertically downward. Probe entrance is 
1 meter from engine accelerator grid. Maximum swing of probe 
is + 18 degrees. 



Figure 8. 


Pulsing Arrangement for the TOF Probe. 


Trigger pulses are generated 60 times per second. Measure- 
ments are made during the first several microseconds following 
the positive step in deflection plate voltage. 
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t = 0 



CAPACITIVELY COUPLED 
TRANSIENT FROM 
DEFLECTOR PLATE 
OFF PULSE 


Hg + ARRIVAL 
~ 2 y SEC RISE TIME 


Figure 9a. Ion TOF Collector Signal. Oscilloscope 
Sensitivity 0.020 Volt/Division Across 
1 Kfl Vertical and 2 ySec/Division Horizontal. 



Figure 9b. Ion TOF Collector Signal With No Ion Beam. 

Capacitively Coupled Transient Still Evident 
But Signal Noise at Later t (>10 Sec) is 
Less Than Trace Width. 



Figure 9c. Ion TOF Collector Signal. Oscilloscope 

Sensitivity 0.010 Volt/Division (Trace Group A) 
and 0.001 Volt/Division (Trace Group B) Across 
In Vertical and 2 pSec/Division Horizontal. 
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IV. Hg ++ /Hg + EXPERIMENTS 

A. ++/+ RATIO AS A FUNCTION OF ANODE VOLTAGE AND ION BEAM CURRENT 

From the previous qualitative analysis of Section II, the most 
important parameters in the ++/+ ratio may be expected to be the anode 
voltage (which establishes the nature of the discharge electron energy 
distribution) and the electron density in the discharge. Measurements 
of the ++/+ ratio as a function of anode voltage are comparatively 
simple to obtain. In contrast, direct measurements of discharge electron 
density are difficult to obtain and assess since they require insertion 
of probes within the discharge chamber with resulting alterations of 
the discharge properties. Therefore, alternative means of establishing 
++/+ dependence on discharge electron density was sought. The method 
which appeared most promising was to examine the ++/+ ratio as a function 
of ion beam current. 

The ion beam current is roughly given by 


■b 4 ‘"Vd 5 td A sg 


( 25 ) 


where I b is ion beam current, qn +d is ion density In the electron bom- 
bardment discharge, v +d Is effective ion drift speed toward the screen 
grid and A g g is screen grid aperture area. Since the plasma discharge 
is highly neutralized throughout most of its volume (n + = n e ), then I b 
is proportionally related to n g . 

Furthermore if I b may be varied without varying the anode potential 
(for example, through variations in propellant injection rate), then 
variations in I b presumably occur without significant shifts in effective 
ion drift velocity, v +d , or in the acceptance area or ion transmission by 
the screen grid. Hence 


I K = k n ,, 

,J v anode fixed 


( 26 ) 


provided that the restrictions outlined above are upheld. 
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Figure 10 illustrates the ratio of ++/+ ion currents in the ion 
beam as beam current and anode voltage are varied. Figure 11 is a cross 
plot of the data in Figure 10. Since the current contribution of an Hg ++ 
ion is twice that of an Hg + ion, the ++/+ current ratio is a factor of 
two larger than the ++/+ abundance ratios utilized throughout the analyses 
[Equations (12) through (18)]. Conversion between the two ratios is con- 
sidered sufficiently simple so that both current and abundance expressions 
of these ratios will be retained. 

Figure 10 has several prominant features which merit discussion. 

The first of these is the "appearance potential," in anode voltage, for 
Hg + production. The double ion content is clearly evident at anode 
voltages of 25 volts and from extrapolations in the cross-plotted data 
of Figure 11, Hg ++ becomes measurable at anode potentials of ^21 to 23 
volts. Since it is unlikely that electrons in the discharge have kinetic 
energies in excess of the cathode to anode potential difference, the 
presence of Hg at these anode settings clearly evidences sequential 
ionization as a contributor to the overall Hg ++ production. (The thres- 
hold for e + Hg° — Hg ++ + e + e is ^30 electron volts, while 
t + Hg + — Hg ++ + e + e has a threshold of ^19 electron volts.) Pawlik 
et al 4 also observed a threshold anode potential of ^21 volts for Hg ++ , 
and concluded that ionization to Hg ++ occurred from other than the Hg° 
ground state. This possibility is discussed below. The data of Milder 1 
also evidenced an onset of Hg ++ at anode voltages of -v20 volts. 

Two possible two-step reactions for the production of Hg ++ may be 
advanced. The first of these is sequential ionization (e + Hg 0 — Hg + + e + e, 
e + Hg + — Hg ++ + e + e). It is also possible to produce Hg ++ through 
ionization of metastable Hg°. THe 6 J P 2 and 6 J Pg states are metastable 
with excitation levels of 5.4 ev and 4.7 ev, respectively, which may be 
subsequently doubly ionized by a single electron impact. The energy 
threshold for this reaction is ^25 electron volts. However, the magnitude 
of processes involving removal of two electrons by a single electron 
Impact remains very small until electron energy Is considerably above 
threshold (see, 11 for example, e + Hg°— Hg ++ + e + e), so that Hg ++ from 
a single stage charge transfer from metastable Hg° is unlikely at 25 volts, 
and, energetically Impossible (for the 6 P states) at 21 volts. In surrenary 
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Figure 11. Doubly to Singly Charged Ion Current Density Ratio as a 

Function of Anode Potential with Beam Current as a Parameter 


to this point, the threshold data in Figures 10 and. 11 are strong 
indications for the production of Hg ++ by a sequential ionization 
process from Hg° to Hg ++ . 

A second feature to note in Figure 10 is the behavior in J ++ /J + 
as a function of beam current for fixed anode voltage. From Equation 
(26), I b is proportional to n g under certain discharge restrictions. 

Also, from Equation (14), the abundance of ++/+ is essentially linear 
in n g provided that n g F^ 2 T is ttie predominant term in the numerator. 

Thus, the ++/+ ratio might be expected to be essentially proportional 
to I b . The data in Figure 10 clearly reveal abundance ratios with, 
essentially, linear dependence upon 1^. As pointed out in the discussion 
of Equation (14), experiments revealing a ++/+ dependence on n g at fixed 
electron energy are strong evidence for sequential ionization. 

A final point to examine in the data in Figures 10 and 11 is the 
dependence of the ++/+ ratio, as anode potential is varied for a fixed 
level of beam current. As may be noted, the ++/+ ratio increases with 
increases in anode potential for all levels of beam current. For fixed 
beam current, the values of n g and T in Equation (19) may be expected 
to remain relatively invariant, yielding the further approximation 


3 anode 


I b fixed 


k I 


( 3F 12 




/ I b fixed 


(27) 


where k is a constant and I b is beam current.* In examining the data 
of Figure 11, the J ++ /J + ratio would appear to indicate that 

^12 

jry -22 volt's") approximately a constant so that one may write 

'anode ' 


< ♦+/♦> =' k ‘ I b (V anode ‘ 22 > < 28 > 

where again, k * , is a constant. It should also be emphasized that the 
analysis here has focused entirely on F-^ as the process through which 


*Note, n g W « 1 . 


28 



the doubly charged ion appears. This is considered a valid assumption 
in view of the estimates of Fq£ derived from the (limit n g — * 0) process 
previously described. 

B. ++/+ RATIO AS A FUNCTION OF ION ACCEL-DECEL VOLTAGE 

The experimental data points given in Figures 10 and 11 were 
obtained with a voltage of +2 kv on the screen grid and with -1 kv on 
the accelerator grid. The plasma thrust beam potential was 0 volt so 
that accelerated ions possessed 2 kev of kinetic energy. In principle 
the accel and decel voltages should not be of consequence in the ++/+ 
ratio. Ions formed in the discharge drift to the forward sheath of the 
plasma and are then accelerated by the applied voltage. Acceleration 
fields do not penetrate into the discharge plasma and, hence, should not 
affect discharge processes. The extent of ion deceleration following 
acceleration should also be of no consequence to discharge processes 
since the decel fields also do not penetrate into the electron bombard- 
ment plasma. To test this concept, the J ++ /J + current ratio was 
determined for a constant 3 kv potential between the screen grid and 
the accel grid as the accel grid potential was varied from -1.0 to 
-2.0 kv, thus varying the ratio of final ion energy to decel potential 
from 2.0 to 0.5, a variation of a factor of 4. The results of this 
experiment are given In Figure 12. The measured current density ratio 
remained essentially constant at 0.12 as expected. Moreover, it is in 
agreement with the results of Pawlik et al 4 for similar variations In 

V and V In summary to this experiment, abundances of ++ and + ions 
" > # 

are not dependent upon ion deceleration fields. 

C. ++/+ RATIO AS A FUNCTION OF DIVERGENCE ANGLE FROM THE THRUST BEAM AXIS 

The experimental data points of Figures 10, 11, and 12 were all 
obtained with the ion T0F analyzer on the axis of the thrust beam. Through 
its mounting arms the T0F analyzer may be moved through the beam. During 
this motion, the T0F analyzer axis continues to be directed at the center 
of the ion engine exit plane so that the entering ion flow remains directed 
along the analyzer axis. The test chamber limits the motion to angles of 
+ 18 degrees with respect to the beam axis, but this incorporates a very 
large fraction of the total ion beam. When J ++ /J + was measured as the 
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analyzer moved from the beam axis to angles of 15 degrees from the beam 
axis, the current ratio remained invariant (Figure 13). This result was 
an expected one, since the emergence and distribution of doubly charged 
Ions relative to the beam axis should be identical to that of singly 
charged ions, as discussed below. 
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If one examines the trajectory of both charge species in moving from 
the plasma discharge to the neutralized beam flow, the radius of curvature 
of a specific charge species is given by 


BEAM DISPERSION ANGLE, 0 (DEGREES) 

Figure 13. Doubly to Singly Charged Ion Current Density Ratio as a 
Function of TOF Analyzer Position Angle 





where q is particle charge, is the potential in -the discharge at 
which the ion originated, V is the potential at the point in question, 
and E x is the electric field perpendicular to the trajectory at the 
point in question. Since q appears in both the numerator and denominator 
of Equation (29), the particle trajectory is charge invariant. This 
trajectory invariance has been previously treated in Section III dealing 
with transverse field deflection in the ion TOF analyzer. Thus, measure- 
ments made on the axis of the beam are representative of the ion beam as 
a whole, as demonstrated in Figure 13. 

D. ++/+ RATIO AS ELECTRON BOMBARDMENT DISCHARGE PARAMETERS ARE VARIED 

In Section IV-A experimental results were presented and compared to 
the predictions of a qualitative analytical model of the discharge. In 
general, agreement is comparatively good, and suggests that at least a 
parametric description of the discharge operation may be obtained from 
the simple analytical model. The discussion has also presented some 
cautionary emphasis regarding use of the simple analysis. This is 
particularly true if a specific variation in the discharge should affect 
the population of more energetic electrons in the discharge. Consider- 
ing that the principal source of doubly charged ions is through sequential 
ionization, processes that affect F-j 2 must be treated with some care. 
Recalling Equation (8) 

F 12 ■ J[\W v e 5 12< £ ) dE < 8 > 

it may be seen that F 12 is sensitive to the "overlap" between f and 
a- 1 2 and, since transitions from the first to the second ionized state 
require electrons with energies in excess of 19 electron volts, discharge 
processes that prevent significant electron populations above 19 ev may 

drastically alter the production of Hg . Such variations that affect 

++ + 

Hg production may not be significant in Hg production since overlap 

between f and Oq-j requires only an electron population above , vlt v ev. 

To illustrate the variations that may be introduced into the ++/+ 
ratio, a series of three experiments was performed with the 20-cm-diameter 
thruster. For all three conditions the screen voltage was maintained at 
+2000 volts with the accelerator grid at -1000 volts. Mercury throughput 
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to the engine was held essentially constant, and since a high value of 
propellant utilization yields the beam current of 0.75 ampere in Case 1, 
below, the utilization in Cases 2 and 3 is substantially reduced. 


Table 1. Variation of Propellant Utilization 


Case and Conditions 


'beamf^P) J V J + 


1. /V 
I 




2 . 



anode 

anode 

magnet 

'anode 

^anode 

^magnet 

anode 

anode 

magnet 


35 volts 

7.5 amperes 

= 0.65 ampere/ 

34 volts 

2.5 amperes 
=0.65 ampere/ 

35 volts 
5.2 amperes 
=0.25 ampere/ 


0.75 


0.45 


0.55 


0.1780 0.0163 >0.90 


0.0047 0.0013 


»0. 55 


0.0138 0.0032 «0.70 


Column E in Table 1 is the standard variation for the experimental 
values of J ++ /J + listed for each case. From the results obtained it 
is evident that large variations in J ++ /J + may occur, under certain 
circumstances, for only modest variations in thruster operating 
conditions. 

In the discussion of a previous section, (IV-A), an approximate 
relationship of 


< ++ / + > 8 k ' >b < v anode ' 22 > (28) 

was introduced as representative of much of the experimental data. To 
illustrate the adequacy of this expression, values of 

^ “ ^++/ J +)/^beam ^anode “ 2 ^) 

are given in Figure 14 for beam currents ranging from 0.25 to 0.90 ampere 
and anode voltages from 25 volts to 55 volts. It may be noted that varia- 
tions in k' in Equation (28) are comparatively small while the quantity 
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Ib (V anoc j e - 22) varies well over one order of magnitude. Even those ob- 
served variations in the calculated value of k 1 could be reduced by varia- 
tions in the magnitude of the chosen "threshold anode voltage" (the 22 
volts used above) by only one volt or so. However, the figure of 22 volts 
drawn from apparent thresholds in Figure 11) is not defined within 1 or 2 
volts. Equation (28) is only a parametric expression of observed engine 
behavior, and minimizing variations in some apparent k' in Equation (28) 
is beyond the accuracy or intent of the analysis. What is of importance 
is to recognize the generally valid linear relationship of J ++ /J + to the 
product of beam current times an anode voltage minus a threshold potential. 

When the data points of Cases 1,2, and 3 of Table 1 are placed on 
Figure 14, it may be noted that Case 1 is in general agreement with the 
"conventional" relationship of J ++ /J + and Ib( v ano d e “ 22 ^ whi ^ e Cases 2 
and 3 are markedly different. It will be advanced here that these 
variations are the result of variations in f under conditions of 
essentially fixed anode potential. 

To produce the operating discharge of Case 2, all engine parameters 
remained at fixed values except that anode current (electron bombardment 
current) was reduced to 2.5 amperes. Nominally this discharge current is 
an order of magnitude larger than beam current, a condition which insures 
a "thick" discharge to single ionization and a high propellant utilization. 
When this anode current diminished to values of only ^5 times beam current, 
the propellant utilization was "spoiled" (diminished to ^0.5). However, 
while 1^ ( v anoc ie ‘ 22 ) only diminishes by a factor of 2 from Case 1 to 
Case 2, the J ++ /J + ratio diminished by ^40. The only apparent source 
for such radical behavior would appear to be that the population of f 0 
above the threshold for sequential ionization has been altered from 
Case 1 to Case 2. While n g and f can remain sufficiently high so that 
approximately half of the injected neutrals are moved from the ground 
state to the first ionized state, the overlap between f and a 12 can vary 
dramatically because of the major variations in a- | 2 moving away from 
threshold. 
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TABLE 1 


Value of Proportionality Constant k 1 in J++/0+ = k ' Ib( V A "22) 
Values are for I^ ranging from 0.25 to 0.90 ampere and V/\ 
ranging from 25 to 55 volts, k' is approximately 0.015 + 

0. 002 for Ib(VA - 22) ranging from ^2 to ^24 watts. Cases 

1, II, III represent other discharge conditions, both 
"on-line" (I, high propellant utilization) and off-line 
(II and III, spoiled propellant utilization) operational 
conditions. 



In moving from Case 2 to Case 3, the anode current was increased 
from 2.5 to 5.5 amperes and magnet current was diminished from 0.65 
ampere to 0.25 ampere. These effects might be expected to be counter- 
acting in their effects on the population of f in energetic electrons, 
since increases in electron bombardment current should increase the 
abundance of energetic electrons in the discharge while diminutions in 
magnetic fields should allow a more rapid transport of fast electrons 
to the anode with a subsequent loss of population of f in energetic 
states. If such counteraction is, indeed, present, it may explain the 
observed values of doubly charged ions, since J ++ /J + remains an order 
of magnitude below the "conventional" level, even for relatively high 

values of discharge parameters < •anode = 5 ' 5 ampcres > V anode ' 22 = 13 
volts) . 

While the results of Cases 2 and 3 and the explanations advanced 
here do not explore exhaustively or explain rigorously all discharge 
behavior, they do suggest that the value of the ++/+ ratio may be 
utilized as a sensitive "probe" of the overall discharge. In this 
Dossible approach to discharge diagnosis, the injected neutrals act as 
probes, with the final abundance ratios of singly and doubly charged 
ions being a complicated but sensitive integral of the "history" of the 
"probing" particle in its passage through the discharge. Possible 
extensions of such techniques could utilize additional "tracer" elements, 
specifically selected for their properties for onsets at various Ioniza- 
tion levels and electron energies. While the present program has not 
pursued these possibilities in depth, it has attempted one further 
experiment in J ++ / J + as anode current is varied for fixed anode voltage 
and beam current. 

Figure 15 illustrates d ++ /J + as I anoc j e decreased from 7.5 amperes to 
3.0 amperes, and with V anQde at 40 volts. Through the region 4.5 amperes 
£ 1^ £7.5 amperes, J ++ /J + remained fixed at M).15 [k ' of Equation (28) = 
0.017 watt"^ in general agreement with the "on-line" operation of Figure 
14]. However, as I anoc j e diminishes to values less than MO 1^, the 
discharge appears to "cool" sufficiently so that ++ production is 
inhibited. This observed behavior is also consistent with Case 2 of 
Table 1 (where the V A = 34 volts is somewhat lower than the V A = 40 volts 
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Figure 15. The J++/J+ Ratio as a Function of Discharge Current. 

Beam current, anode potential, and other operating 
parameters were held fixed. 

of Figure 15). From Figure 5, it is apparent that good propellant 
utilization may be obtained. with diminishing fractions of doubly charged 
ions, by suitable reductions in anode current. Below certain discharge 
current levels, however, the discharge becomes noisy and may be unsatis- 
factory for other reasons relating to ion engine impact upon spacecraft 
performance. 

The data of Figure 5, and that of Figure 10 and Table 1 are also of 
interest to the point that they suggest that electron density in the ion 
engine discharge is not closely coupled to electron bombardment current 
(at least, for high propellant utilization discharge conditions). If 
beam current is, indeed, related to discharge electron density as sug- 
gested in Equations (25) and (26), then electron density in the discharge 
for the data of Figure 15 remained essentially constant as anode current 
varied over a factor of ^2.5 (here propellant injection remained fixed). 
In Figure 10, variations in beam current of a factor of 4 were obtained 
with anode current and voltage fixed and variations in the neutral 
injection rate. 
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A final area for comment on the data in this and the preceding 
portions of Section IV is in the apparent threshold for ++ production. 
From Figures 10 and 11, it would appear that anode voltages of ^22 volts 
are required for the onset of ++ production. Presumably electrons in 
the discharge may acquire kinetic energies as large as e V anoc j e > However, 
electrons in the discharge are first present in the region between the 
hollow cathode and keeper and then move into the main discharge chamber. 
In a separate series of experiments investigating electron temperatures 
in the plumes of hollow cathode plasma bridge neutralizers, it has been 
found that electrons have energies which are small compared to e(V|< ee p er 
V ca thode^' ^h 1s su 99 es1:s that the P^ asma densities in these hollow 
cathode discharges are sufficiently dense that electrons are in collision 
ally dominated transport, losing energy to collisions with atoms and ions 
and gaining energy from the electric field. In this manner an electron 
leaving a surface at V cathode eventually dwell in a plume at V keeper 
and may only acquire a kinetic energy which is small compared to the 
indicated differences in potential. In this same manner, electrons in 
the ion engine discharge may reach the keeper region with only thermal 
energies of fractions of ev. Their subsequent entry into the main dis- 
charge may be sufficiently collisionless to allow rapid acquisition 
of kinetic energy but the maximum possible electron energy would then 

bease(V anode " V keeper ) ’ When the keeper P otentia1 of 3 volts is 
included, then the earlier expression in Equation (28) utilizing 

(V anode ‘ 22 > v0,ts ™ y be rewr1tten as [< v a node ' v keeper> * 19 3 volts ’ 
which may be again rewritten to C(V anode . V keeper ) - V 12>threshold ] 

where V 12 threshold 55519 volts is the known energetic threshold for 

sequential ionization of Hg + . It has been previously suggested by 
1 ? 

Pawlik that the important parameter in the study of ++ production may 

be < V anode " Vkeeper 1 * The apparent onset of ++ production for V anode " 

V keeper “ V 12, threshold P resents interesting possibilities, both in 
explanation of observed ++ production, and, (referring to previous 
considerations of ++ production as a "probe" of plasma discharge 
behavior) as an element of ion engine discharge diagnosis. 
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V. DISCUSSION AND SUMMARY 


The ion TOF probe results have indicated a threshold anode potential 
of approximately 22 volts for the appearance of Hg ++ . This is in quali- 
tative agreement with the early data of Milder^ for a somewhat different 
form of electron bombardment discharge, and is in good agreement with 
the data of Pawlik et al^ for an identical 20-cm bombardment thruster. 

The 22-volt appearance potential for Hg ++ is strong evidence that the 
doubly charged ion results from sequential ionization, Hg°— Hg + — Hg ++ . 

The magnitude of Hg ++ /Hg + observed in these experiments is in good 
agreement with the results of Pawlik et al for those ions emerging near 
the axis of the ion thruster. Relationship of the response of the two 
probes is given in Appendix A, which follows. From the treatment there 
it is shown that the ion TOF probe results should be compared to Pawlik's 
on-axis data. 

The Hg ++ /Hg + ratio has been shown to depend in an approximately 

linear fashion on n . the electron density in the discharge. This, like 

6 

the threshold data, is strong evidence for the sequential ionization 
process for Hg ++ production. [Any two-step process, for example 
Hg°— Hg°* (metastable) — Hg ++ , should also possess a linear dependence 
on n . However, the magnitude of the cross section for 0* to 2 + transl- 
tions makes this reaction an unlikely source of Hg ]. 


The ++/+ ratio has been examined and found to be invariant with 
respect to variations in the ion deceleration fields following acceler- 
ation and with respect to variations in ion divergence angle relative to 
the thrust beam axis. Both aspects of this behavior are expected. 


The ++/+ ratio has been found to be parametrically described by 
k I b ( V a " ^2) where 1^ is ion beam current and (V^ -22) is anode voltage 
minus 22 volts and is approximately [(V A - V keeper ) - v thre|hold^ where 
^threshold t ^ ie a PP earance potential for ionization of Hg to Hg ’ 

This relationship has been observed to hold for high propellant utilization 
conditions. However, in a "spoiled" discharge with poor propellant utili- 
zation conditions, marked variations of the ++/+ ratio away from this 
stated parametric dependence have been observed. This behavior has also 
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been examined as anode current was varied. Within a range of anode 
current from approximately seven to approximately fourteen times ion 
beam current, the ++/+ ratio remained fixed. 

The apparent decoupling between anode current and discharge plasma 
density (variations of 1^ through variations in propellant flow and with 
fixed anode voltage and current, and constancy of I b with varying anode 
current) raises interesting questions. Kaufman and Cohen^ state that 
primary electron density is linearly dependent on anode current. However, 
to compare results from these experiments to Kaufman and Cohen may require 
a fixed expenditure in ev/ion as other parameters in the discharge are 
varied. For those experiments in which 1^ was varied through neutral 
injection rate, the ev/ion cost remained above the 300 to 400 ev/ion 
figure but did vary inversely with I b - Similarly, for fixed I b and 
varying anode current, the ev/ion cost remained above 300 to 400 ev/ion 
point, but varied directly with anode current. 

A final aspect of these experiments is the possibility they reveal 
of using the doubly charged ion abundance as a probe of plasma discharge 
behavior. Admittedly, the ++/+ ratio is an involved result of the ion 
history as the particle moves through the discharge, and does not provide 
a point-by-point diagnosis of that discharge. However, the ++/+ ratios 
are easily and quickly obtained and do provide an indication of discharge 
behavior without the (possibly disruptive) insertion of material probes 
into that discharge. 
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APPENDIX A 


ESTIMATE OF PROBE SENSITIVITY FUNCTION 

In comparing the results of Hg ++ /Hg + ratios obtained here to those 

4 

obtained by Pawlik et al , attention must be given to acceptance charac- 
teristics of the two probes. The Pawlik probe accepts a well -collimated 
ion beam and, hence, "looks" at only a small portion of the ion engine 
accel surface at a given probe position. Results of Pawlik, then, 
describe the ++/+ ratio as a function of radial exit position on the 
accel grid, and a diminution of this ratio as the probe viewing point 

moves from the axis of the thruster to the outer edge of the discharge 

• 4 

region has been observed. 

The transverse field ion TOF probe has a larger acceptance range 
on the ion entrance angle, and, qualitatively speaking, "looks" at the 
entire "active" face of the accel grid. However, certain portions of 
the accel grid are viewed with more sensitivity than others, because of 
geometrical factors in the probe inlet channel and sensing Faraday cup. 

The probe entrance aperture is a 1- by 3-cm rectangular opening followed 
by a series of deflecting plates extending about 4.5 cm into the probe 
housing. The accel grid-to-probe entrance distance Is ^1 meter and the 
probe entrance-to-Faraday cup separation is ^0 cm. Ions originating on 
the axis of the Ion engine and passing through the entrance aperture will 

illuminate an area on the Faraday cup of (1 x 3) cm. 

This area is totally contained within the limits of the Faraday cup, 
and the inlet solid angle acceptance is at maximum value so the "sensitivity" 
to ion detection is evaluated at 1.00. Within a narrow range of area over 
the face of the accel grid this sensitivity remains at 1.00. Outside of 
this area, ions have a diminished solid angle of acceptance through the 
probe inlet. Also, the illuminated area at the plane of the Faraday cup 
may not be completely contained within the sensing area of the cup. When 
these two geometrical factors are evaluated, a sensitivity function for 
the entire active area of the accel grid Is determined. This function is 
given in Figure 16. 
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Figure 16. A Ray-Tracing Estimate of the Sensitivity Function of the 
Time-of-Fl ight Probe to Various Regions of the Ion Engine 
Grid Normalized to Center 

The variation in sensitivity for ion acceptance over the face of 
the accel grid will not result in variations in the measured ++/+ ratio 
if the ++/+ production ratio in the discharge is invariant over this 
same area. However, if the ++/+ ratio should diminish in moving away 
from the engine axis, and Pawlik's results indicate such a diminution, 
then corrections of the experimental data from the TOF probe must be 
made. Such corrections involve the calculated sensitivity function and 
the measured dependence of J + over the accel grid. Expressed in 
integral form: 



(A-l) 
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S dA 


+ / 


measured 


f (±) 

beam area'^' 


(A-2) 


S dA 


where — J , J , and S (the sensitivity function), are all functions 

of position over the area of the beam. In the absence of the measured 
/J + \ 

behavior of ( — ) over the beam area, the correction from Equation (A-2) 

to Equation (A-l ) cannot be made. It should be noted, however, that S is 

/ J +\ 

generally large where I — J has its maximum values for this bombardment 

discharge, and, hence, the measured ++/+ ratios are accurate representations 
of the overall beam specie currents and of the ++/+ ratios on and near the 
axis of the ion engine. In making comparison with Pawlik's results, then, 
the TOF probe will be compared to the on-axis results of that probe 
(Figure 11, Reference 4). 
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